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We use the decay modes B — > .Ko(1430)7r and B — ► ao(980)ii' to study the scalar mesons Kq (1430) 
and ao(980) within perturbative QCD framework. For B — » Kq (1430)-7r, we perform our calculation 
in two scenarios of the scalar meson spectrum. The results indicate that scenario II is more favored 
by experimental data than scenario I. The important contribution from annihilation diagrams can 
enhance the branching ratios about 50% in scenario I, and about 30% in scenario II. The direct CP 
asymmetries in B — * ifo(1430)7r are small, which are consistent with the present experiments. The 
predicted branching ratio of B — > ao(980)A" in scenario I differs from the experiments by a factor 2, 
which indicates ao(980) can not be interpreted as qq. 



I. INTRODUCTION 

The scalar meson spectrum is an interesting topic for both experimental and theoretical studies, but the underlying 
structure of the light scalar mesons is still under controversy. Many scalar meson states have been found in experiments: 
isoscalar states cr(600), /o(980), /o(1500), /o(1370), /o(1710); the isovector states ao(980), ao(1450) and isodoublets 
k(800), A*q(1430). In the literature, there are many schemes for the classification of these states [1-4]. Here are two 
typical scenarios: the members of the lower mass nonet cr(600), k(800), /o(980), and ao(980) are treated as the lowest 
lying qq states, while ATo(1430) et al. which form the higher mass nonet are the first excited qq states; In scenario 
II, the members of the lower mass nonet are treated as the four-quark states [1]. Then the higher mass nonet is 
considered as the lowest lying qq states. There are also other schemes to classify these states, for example, c(600) 
and k(800) are not considered as the physical states, a (980) (or a (1450)), / (980), A"q(1430) and / (1500) form the 
qq nonet [3]. In this paper, we study the scalar mesons in the first two scenarios. 

Although intensive study has been given to the decay property of the scalar mesons, the production of these mesons 
can provide a different unique insight to the mysterious structure of these mesons. Compared with D meson decays, 
the phase space in B decays is larger, thus B decays can provide a better place to study the scalar resonances. 
Experimentally, B meson decay channels with a final state scalar meson have been measured in B factories for 
several years [5]. Much more measurements have been reported by BaBar and Belle [6-10] recently (see [11] for more 
experimental data). On the theoretical side, the B decays which involve a scalar meson have been systematically 
studied using QCD factorization (QCDF) approach by Cheng, Chua and Yang [12]. They draw the conclusion that 
scenario II is more preferable. For example, in scenario I, the predicted branching ratio of B~ — > Kq°it~ without 
annihilation is only about 1 x 10~ 5 , which is much smaller than the experimental results. In order to explain the 
large data, the annihilation contribution is required to be large. But in this scenario, a (980) is also a qq state 
and the SU(3) symmetry implies a much large annihilation contribution to B — > ao(980)K . This large annihilation 
contribution can lead to a larger branching ratio than the experimental upper bound. Then it is concluded that 
scenario I is less preferable than scenario II. 

The annihilation topology contribution plays such an important role that we should pay much more attention to it. 
In QCD factorization approach, which is based on collinear factorization, the so-called endpoint singularity makes the 
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annihilation contribution divergent. This might be resolved by applying the "zero-bin" subtractions which will lead to 
new factorization theorems in rapidity space [13]. The more popular way to handle this divergence is the Perturbative 
QCD (PQCD) approach [14, 15]. Using this approach, some pure annihilation type decays have been studied and 
the results are consistent with the experiments [16], which indicates it a reliable method to deal with annihilation 
diagrams. In the present paper, we will use the PQCD approach to calculate the decay modes B — > ^(1430)^ 
and B — > ao(980)i4r (in the following, we use Kq and ao to denote JTq ( 1430) and ao(980) for convenience). We will 
examine how large the annihilation topology contribution is, and find whether there are enough reasons to determine 
which scenario is more appropriate through these decay channels. 

This paper is organized as follows: in the next section, we give a brief review of the scalar meson wave functions, 
which are important inputs in PQCD approach. What followed is the analysis of the B — > KqTt and B — ► a K decays. 
The numerical results and the discussions are given in the section IV. Our conclusions are presented in the final part. 



II. SCALAR MESON DISTRIBUTION AMPLITUDES 

In the two-quark picture, the decay constants fs and fs for a scalar meson S are defined by: 

(S(p)|927^i|0) = fsP^ (S(p)|?2<7i|0> = fsms, (1) 

where ms(p) is the mass (momentum) of the scalar meson. The vector current decay constant fs is zero for neutral 
scalar mesons a, fo and a° due to the charge conjugation invariance or G parity conservation. In the SU(3) limit, the 
vector decay constant of K$ also vanishes. After including the SU (3) symmetry breaking, it only gets a very small 
value which is proportional to the mass difference of the constituent quarks. The scalar density decay constant fs 
can be related to the vector one by the equation of motion: 

Ms/s = fs, (2) 

where us is defined by ris — r, ul^L — rrr? which is scale dependent, thus the scalar decay constant is also scale 

dependent. Many model calculations have been performed for the Kq and ao [17]. In this paper we will use the values 
from QCD sum rules in ref. [12]. Fixing the scale at lGeV, we specify them below: 

scenario I : Ik* = -(300 ± 30)MeV, fx* = -(25 ± 2)MeV, 

f ao = (365±20)MeV; (3) 
scenario II : f K * = (445 ± 50)MeV, f K * = (37 ± 4)MeV. (4) 

Now we turn to the distribution amplitudes. The scalar meson's light-cone distribution amplitude is defined by: 
(S(p)|5«(z)«2„(0)|0) = i / dxe"'-Ul>s(*) + ms4(n) - ^svWSW 



±f l 



a/3 

dxe ix P- z {^cl>s(x)+ms(t> S s(x)+msm-l)^(x)} , (5) 

a/3 



where n = (1, 0, 0t) and n — (0, 1, 0t) are dimensionless vectors on the light cone, and n is parallel with the moving 
direction of the scalar meson. The distribution amplitudes 4>s {x), <Ps( x ) an( i 4>s{ x ) are normalized as: 

J^dx^x) = J\x4(x) = J^dx^ix) = A, (6) 

and <%{x) = \£ i <tr s (x). 
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In general, the twist- 2 light cone distribution amplitude 4>s{%) can be expanded as: 

oo 

Mm) 



<t>s(x,(J.) = ^=-6x(l-x) 



2\/6 



■6x(l — x) 



m—l 

00 



l + ^ s J2B m ( (J ,)C 3 J 2 (2x-l) 



m—l 



(7) 



where B m (p) and Cm 2 (a;) are the Gegenbauer moments and Gegenbauer polynomials respectively. The Gegenbauer 
moments B\, B3 of distribution amplitudes for Kq and ao have been calculated in [12] as 

scenariol: B x = 0.58 ± 0.07, B 3 = -1.20 ± 0.08; (8) 
B x = -0.93 ±0.10, B 3 = 0.14 ±0.08; (9) 
scenarioll: B l = -0.57 ± 0.13, B 3 = -0.42 ± 0.22, (10) 

where the second line is for a , and the others are for Kq. These values are also taken at fi = lGeV. 

As for the twist-3 light-cone distribution amplitudes, there is no study on their explicit Gegenbauer moments so 
far, so we take the asymptotic form in our numerical calculation: 

= <i> T s(x) = ^L(l-2x). (11) 

In our calculation, we will choose the momentum fraction on the anti-quark, thus we should use </>s(l— x), (j)g(l—x) = 
4> s (x) and </>g(l — x) = — 4>^{x). But in the amplitudes, for simplicity, we use 4>s(x) to denote </>s(l — x) . It is similar 
for the pseudoscalar meson. 

III. THE PERTURB ATIVE QCD CALCULATION 

In this section we will give the decay amplitudes for the B — > Kq-k and B — > a K decays in the PQCD approach. 
The PQCD approach is based on the kr factorization [18], where we keep the transverse momentum of the partons 
in a meson. In this approach there is no divergence when the longitudinal parton momentum fraction falls into the 
endpoint region. The weak decay matrix element can be completely factorized [19]: 

A = <p B ® H (6) ® J (g) S (g> (j> Ml ®4>m 2 , (12) 

where (f>B, (j>M 1 and (j)M 2 denote the wave functions of the B meson and the light mesons, respectively. S and J denote 
the Sudakov form factor and the jet function respectively. The Sudakov form factor comes from kr resummation which 
kills the end-point singularities. The jet function is from threshold resummation which can organize the large double 
logarithms in the hard kernel. The symbol £g> denotes convolution of the parton longitudinal momentum fractions 
and the transverse momentum. is the six-quark hard scattering kernel, which consists of the effective four quark 
operators and a hard gluon to connect the spectator quark in the decay. The standard four-quark operators describing 
the b — > s transition are defined as [20] : 

• current-current tree operators 

Oi = (u a bp)v-A(Sl3U a )v-A, 2 = (u a b a ) V -A(Sl3U0)v-A, (13) 

• QCD penguin operators 

Os = {s a b a ) V -A^%q'p)v-A, O A = {s b a ) V - A ^{ti a q'p)v-A, (14) 

q> q> 

5 = {s a b a ) V -A^{q'pq'p)v+A, O e = (s(ib a )v-A^2(q' a q'f3)v+A, (15) 

q> q> 
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• electro-weak penguin operators 

3 
2 



3 ^ 3 

07 = -^{s a b a )v-A^e ql {q fj q' p )v+A, Os = jfe^lv-i^^fM^+A, 
<?' <?' 

q' q' 



(16) 
(17) 



where q' = (u,d,s,c,b). The 10 operators together with their QCD-corrected Wilson coefficients form the effective 
Hamiltonian: 

10 "| 



Gf 
V2 



Vu b v: 



Ci(/i)Oi(/i) + C 2 ^)0 2 (fx) - VtbVt'Y^CiMOifa) . 



(18) 



l — 3 



The partition of the perturbative and non-perturbative region (factorization scale) is quite arbitrary, but the full 
amplitude should be independent of the partition. We usually take the largest virtuality of internal particles as the 



factorization scale, which is of order vWA. The leading order Wilson coefficients will be evolved to this scale. As 
the factorization scale in PQCD approach is smaller than the factorization scale of QCDF, which is about me, a large 
enhancement of the Wilson coefficients occurs, especially for the penguin operators [14, 15]. 



A. B — > K*)Tx decays 

The leading order Feynman diagrams for these decays in PQCD approach are given in Fig. 1. The decay amplitude 
for each diagram can be obtained by contracting the hard scattering kernels and the meson's wave functions. According 
to the power counting in PQCD approach [21], the first two emission diagrams in Fig. 1 give the dominant contribution. 
For the (V — A)(V — A) kind of operators, the decay amplitudes for these two diagrams are given by: 

327T f°° f 

Fg_ n (a) = —m B f K - J dx x dx 3 J M&i b 3 db 3 (j)g(x 1 , bi)^a(t)E e (t) 

x [(1 + x 3 )^(x 3 ) + rv(l - 2x 3 ) (<f£(x 3 ) + <t£(x 3 j)] h e (x u x 3 , b l7 b 3 ) 

+2r^(x 3 )a(t')E e (t')h e (x 3 ,x 1 ,b 3 , &i) j, (19) 

where r„ = m^/mg, is the chiral enhancement scale and a is the corresponding Wilson coefficient. E e (t) is defined 



as 



E e (t) = a s (t)eMSB(t)-S v (t)}. (20) 

For the (V — A)(V + A) kind of operators, the decay amplitudes for these two diagrams are given by: 

647r f 00 f 

F B^i a ) = -— m B r K^.fK' Q J dx x dx 3 J b l db l b 3 db 3 4> B (xi 1 b 1 )la(t)E e {t) 

x [^(x 3 ) + r n x 3 (<t>%(x 3 ) - 4>l(x 3 )) + 2r 7T <f>P (x 2 )] h e {x 1 ,x 3 MM) 

+2r^{x 3 )a{t')E e (t')h e {x 3 ,x 1 ,b 3 ,b 1 )^ , (21) 

with the factorization scales t = maxjy^ims, l/b 3 } and t' — max {y^j"ms, 1/fei, I/63}, rjf* = m K */m B . 
The Sudakov form factors Sg(t) and S n (t) and the hard functions h e and others like h a , h na are given explicitly in 
ref.[14, 15]. 

Comparing Fg^^ and Fg^ n , we find that the first one is proportional to the small vector decay constant, but the 
latter is proportional to the scalar decay constant which is strongly chiral enhanced, so Fg^ will give the dominant 
contribution. 
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FIG. 1: The leading order Feynman diagrams in PQCD for B — > Kq-k and B — > aoA' 

In ref. [12], the author found that the vertex corrections and hard-spectator-scattering corrections can enhance (24 
sizably. In PQCD approach, the vertex corrections are at the next-to-leading order in a s , so we neglect it in our 
calculation, but we include the hard spectator scattering (the last two diagrams in the first row of Fig. 1). After the 
calculation, the non- factorization decay amplitudes for the (V — A)(V — A) kind of operators read: 



M L B ^{a) = 



1287T 



3 ^m%J dxidx 2 dx 3 J b\db\ b 2 db 2 4>b(%i, h)(f>K* (x 2 ) 

x| - (x 2 - l)(f) v (x 3 ) + r 7I x 3 (4>^(x 3 ) - <t%(x 3 j) a(t)E' e (t)h n (xx,l - x 2 ,x 3 ,h,b 2 ) 
- (x 2 +x 3 )(f)^(x 3 ) -r„x 3 (0f (x 3 ) + 4%(x 3 )) a(t')E / e (t')h n (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

For the (V — A)(V + A) kind of operators, the decay amplitudes for these two diagrams are given by: 

1287r Z" 00 f 

Mg^( a ) = m 4 B r K * J dx\dx 2 dx 3 J bidbi b 2 db 2 (/>b(xi, bi)yi(t)E' e (t) x 

(x 2 ~ 1) ^(x 3 )(d> s K »(x 2 ) + <t>l s (x 2 )) +r n (tf(x 3 ) - <ft(x 3 ))(cf> s K «(x 2 ) + &„(x2)) 

-r^x 3 (4>^(x 3 ) + <j%{x 3 )) ( cj>K* (x 2 ) - 4>k» (x 2 U K(x u 1 - x 2 ,x 3 , h,b 2 ) + 
x 2 <j>i{x 3 ) (<t> S K .{x 2 ) - <f> T K *{x 2 )} + r n x 2 (<fi(x 3 ) - <j>l{x 3 )) (cj> S K *(x 2 ) - <j> T K .{x 2 ) 

+r n x 3 (4>k(x 3 ) + 4>l{x 3 )) \ 4> S k*{x 2 ) + (I>k*( x 2)) a{t')E' e (t')h n (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 



(22) 



(23) 
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where E' e {t) = a s (t) exp[— S B (t) — 52 (t) — 53(t)], the factorization scales are chosen by 



t = maxl^/iissms, -x x - x 2 )x 3 \m Bl 1/h, l/b 2 }, 
t' = maxl^/iissms, V\( x i ~ x 2 )x 3 \m B , 1/bi, l/b 2 }. 



(24) 
(25) 



From the above formulae, we can see that the two hard spectator scattering diagrams contribute constructively 
for M. B _^ while most contributions cancelled for M B ^. So it is expected that the (V — A)(V — A) kind operator 
contribution can give an important contribution as in B — ► f K [22]. However, these diagrams are suppressed 
compared with the factorizable ones for a smaller Wilson coefficient C 3 , thus they cannot play such an important role 
as in QCDF [12]. 

In PQCD approach, the annihilation type diagrams are free of endpoint singularity, so they can be calculated 
systematically. The Feynman diagrams are plotted in the second row of Fig. 1. For the first two factorizable 
annihilation diagrams, the decay amplitude formulae are written as: 

327T f°° 

Fa (a) = —m B f B J dx 2 dx 3 J b 2 db 2 b 3 db 3 

x j [(x 3 - 1)44(x 3 )<I>k 5 (x2) - 2r v r K .{x 3 - 2)^ {x 3 )<j> s K , (x 2 ) 

+ 2r 7T r K5 x 3 (f)K,(x 2 )<f>l(x 3 ) a(t)E a (t)h a (x 2 , 1 - x 3 , b 2 ,b 3 ) 
+ [x2(f4(x 3 )4>K S (x2) - 2r 7r r K .jP(x 3 ){{x 2 + 1)4k~{x 2 ) + {x 2 - l)^«(x 2 )) 

xa(t')E a (t')h a (l - x 3 , x 2 ,b 3 , 6 2 ) j, 

for the (V — A)(V — A) kind of operators and 

64tt 



(26) 



TO = -—m%f B 



J dx 2 dx 3 J b 2 db 2 b 3 db 3 ^E a (t)a(i)h a (x 2 , 1 - x 3 ,b 2 , b 3 ) 

r v (X3 - 1)4k*(x 2 ) (<Pv(X3) + 4%{X3)) + ^r K ^{x 3 )4 s K ,{x 2 ) 
2r 7I (t) KS (x 2 )(j)^(x 3 ) + r K *x 2 (j) v (x 3 ) {4> T K *{x 2 ) - </>|.(x 2 )) 

xa(t')E a (t')h a (l - x 3 ,x 2 ,b 3 ,b 2 ) \, 



(27) 



for the (V — A)(V + A) kind of operators, where 

E a (t) = a s (t)exp[-S 2 (t)-S 3 (t)}, (28) 

with t = max{-y/l — x 3 m B , l/b 2 , l/b 3 } and t' = maxjy^mB, l/b 2l V^}- 

For the non- factorizable annihilation diagrams, e.g., the last two diagrams in the second row of Fig. 1, the factor- 
ization formulae read: 

1287r Z* 00 f 

Ma(a) = -^=m B J dxidx 2 dx 3 J b 1 db 1 b 2 db 2 4 B (x 1 ,b 1 )l[-x 2 4^(x 3 )4 K *(x 2 ) 

+r 7V r K ,(j)P{x 3 ) {{x 2 - x 3 + 3)(j) S K ,(x 2 ) + (x 2 + x 3 - l)(j) T K , (x 2 fj 
+r v r K *4l{x 3 ) ((1 -X2- x 3 )<j>K. (x 2 ) + (1 - x 2 + x 3 )4%, (x 2 )j 
xa(t)E' e (t)h na (xi,x 2 ,x 3 ,b 1 ,b2) - [{x 3 - l)4*(x 3 )4> K *{x2) 
+r v r K *<j>v{x 3 ) {(X2 -x 3 + l)4 s K ,(x 2 ) - (x 2 +x 3 - 1)4%. (x 2 )j 



+r v r K *(j%(x 3 ) {{x 2 + x 3 - l)4> s K *{x 2 ) - (1 + x 2 - x 3 )4>k,(x 2 )^ 
xa{t')E' e {t')h' na (x 1 ,X2,x 3 ,b 1 ,b 2 )] , 



and for the (V — A)(V — A) kind of operators 

1287T 



M2(a) = 



— i=m% I dx\dx 2 dx 3 I b\db\ b 2 db 2 4>b{x\, b\) x 
3y6 Jo Jo 

r v {x 3 - l)<]) K *{x 2 ){<j>l{x 3 ) - 4>„{x 3 )) + r K *x 2 (j) n {x 3 ){4) S K .{x 2 ) + 4>^{x 2 )) 

xa(t)E' e (t)h na (x 1 ,x 2 , x 3 , h,b 2 ) 
r v {x 3 - l)(j) K 4x 2 ){(j)l{x 3 ) - 4>^{x 3 )) + r K *x 2 <j> v (x 3 ){<t> s K Jx 2 ) + <t> T K Ax 2 )) 



xa(t )E e {t )h na (xi,x 2 ,x 3 ,h,b 2 ) } , 
for the (V — A)(V + A) kind of operators, where 



t = niax{i/s 2 (l - x 3 )m,B, \/l - (1 - xi - x 2 )x 3 itib, 1/&1, 1/&2}, 



t' = maxf^fl - 13)™^ vK^i ^ x 2)(l - x 3 )|m B , 1/62}- 
Summing up all contributions mentioned above, the decay amplitude for B — > K£°Tr~ is 



A(B-^if* l 7r-) 



= Gf 

-V tb V*^F e L (a 4 - loio) + F e fl (ae - ^s) + (04 + do) + F*(a 6 + a s ) 
+M L e (C 3 - l -C Q ) + M^(C 5 - l -C 7 ) + M L a {C 3 - l -C Q ) + M«(C 5 - l -C 7 ) 



where the combinations of Wilson coefficients are defined as usual [23] : 

01 = C 2 + Ci/3, a 3 - C 3 + C4/3, a 5 - C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = C 9 + C w /3, 
a 2 = Ci + C 2 /3, a 4 - C 4 + C3/3, a 6 = C 6 + C 5 /3, a 8 = C 8 + C 7 /3, ai = C w + C 9 /3. 



(29) 



(30) 



(31) 



(32) 
(33) 



For the decays B° — > /fQ _ 7r + , B~ — > ifQ~7r° and _B° — > /QJ tt°, the analysis is similar, except the last two channels 
include the 7r-emission diagrams in the third row of Fig. 1. The decay amplitudes of the factorizable 7r-emission 
diagrams for the (V — A)(V — A) kind of operators are written by: 

32tt f°° ( 

F B^K*i a ) = — m %U J dxxdx 2 J b 1 db 1 b 2 db 2 (j> B {xi,b 1 )^a(t)E e (t) 

x (1 + x 2 )(j) K .(x 2 ) - r K ,(l - 2x 2 ) (<i> s K *{x 2 ) + 4>]i*{x 2 )^ h e (x 1 ,x 2 ,b 1 ,b 2 ) 

-2rK*(j> S K*{x 2 )a{t')E e {t')he{x 2 ,x\,b 2 ,b\) >, 



and for the (V — A)(V + A) kind of operators: 



Ft 



-F. 



For the non-factorizable diagrams (V — A)(V — A) operators: 

1287T 



(34) 



(35) 



1287T f'^ f°° ( 

■M B ^ K *{a) = -^-j=m B J dxidx 2 dx 3 J hdh b 3 db 3 4>b{xu h)<l>n(x 3 )la(t)E' e 



(*) 



(1 - xs)(j)K»{x2) + r K *x 2 \4>k*(x 2 ) - <Pk*(x 2 )J h n (xi,l - x 3 ,x 2 ,bi,b 3 ) 
+ -{x 2 + x 3 )(f> K * (x 2 ) ~ r Ks x 2 ((f> s K , (x 2 ) + <fe, (x 2 )) 
a(t , )E' e (t')h n (x 1 ,x 3 ,x 2 ,b 1 ,b 3 )^ , 
and for (V — A)(V + A) operators 

f 1 f°° ( 

' - J d Xl dx 2 dx 3 J b 1 db 1 b 3 db 3 ^ B (x 1 ,b 1 )^(x 3 )l-a(t)E' e (t) 

h n (xi, 1 - x 3 ,x 2 M,b 3 ) 



(36) 



, , . 1287T 4 

M^ K ,(a) = -to 4 



(x 2 - x 3 + l)<p K * (x 2 ) + r K *x 2 ((/) s K , (x 2 ) + <\> T K , (x 2 )) 



+ 



x 3 4>k*{x 2 ) +r K *x 2 (cp s K ,(x 2 ) - 4>k*(x 2 )) a(t , )E' e (t')h n (x- L ,x 3 ,x 2 ,b 1 ,b 3 )^ . (37) 
We write the decay amplitudes for the other three channels below: 
A(B° K*-n+) =9z[ Vub v: s {F^M+M^iCi)} 

-V tb V*{F^ n (a 4 + a w ) + F*^(a 6 + o 8 ) + F^(a 4 - ±a w ) + F*(a 6 - ^a 8 ) 

+M L B ^(C 3 + C 9 ) + M%^(C B + C 7 ) + M L a (C 3 - \c 9 ) + M*(C 5 - \c 7 )}] , (38) 
A(B- K* -n°) = 9l. [v ub V* s {F^M + M^(d) + F^ K , (a 2 ) + M L B ^ K , (C 2 ) + F a L ( ai ) + M^d)} 

-V tb V*{F^ n (a 4 + oio) + Fg^(a e + o 8 ) + F^ K ,(L 9 - ^a 7 ) 
+F a L (a 4 + oio) + F*(oe + as) 

+M L B ^(C 3 + C 9 ) + M^(C 5 + C 7 ) + M L B ^ K , {\c w ) + M^ K . (|c 8 ) 
+MZ(C 3 + C W ) + M?{C S + C 7 )}], (39) 
A(B° -+ Kfn°) = ^ [V ub v: s {F^ K , (a 2 ) + M L B ^ K , (C 2 )} 

-V tb V t * s { F^Ja 4 ioio) - F^Ja 6 \a s ) + F^ K , (^a 9 ) + (|a r ) (40) 

- ^a 10 ) - F«(a 6 - ^a 8 ) 
-M^(C 3 - ic 9 ) -Mi^(C B - ic 7 ) + A^.(^io) + A<i^.(|c 8 ) 
-MJ(C 3 - ^C 9 ) - M?(C 6 - ic 7 )}] . (41) 
The isospin relation for the four channels holds exactly in these equations: 

V2A(B° -» i^V) + A(B° K*-tt+) = V2A(B- -» i^Tr ) - A(B~ -» K^tt - )- (42) 



B. B — > ao-ff decays 



As mentioned above, the predicted branching fractions of S — > ao-K", which overshoot the experimental limits, are 
regarded as an evidence to rule out scenario I in QCDF, thus it is important to see whether it is the same in the 
PQCD approach. 

The Feynman diagrams for these decays are completely the same as the B —> K^ir except that we should identify 
the sq and qq as K and a 9 rather than Kq and tt, then each channel corresponds to the one in B — > K^tt. Their 
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factorization formulae can be derived from the corresponding channels directly. The if-emission and annihilation 
decay amplitudes of B — > a K can be obtained from B — > A'qTt by making the substitution: 



^ -aS ±P jT ,T 



^^^ ao ,^^-^ o ,^^-^ o . (43) 
The ao-emission diagrams have only nonfactorizable contributions, the substitution for (V — A)(V — A) operators is: 

(44) 



(45) 



but for (V - A)(V + A) operators, 



^->0o„,0£-><,^->C (46) 



(47) 



Compared with f? — > ii"o7r, the features of B — ► aoA' are: 

• For the decays £?° — > an d ^ — * a o tnc emitted particle is if, which can be produced through the 

axial- vector current without any suppression, thus the operator O4 can give a large contribution to the emission 
factorizable amplitudes. But this term has a minus sign relative to Oq, so the penguin operators cancel with 
each other sizably. The contribution from tree operators can be large due to the large Wilson coefficients ai in 

1 



B~ -» a+K~ 



• In B — > ATq 7r~, the emitted particle (K£ ) is a scalar meson, the two hard spectator scattering diagrams 
(non-factorizable) can enhance each other due to the anti-symmetric twist- 2 distribution amplitudes. But in 
B° — ► clqK~ and B~ — » a^K , the emitted particle is a pseudoscalar, there are cancellations between the two 
hard spectator scattering diagrams. So the hard spectator scattering contribution is rather small. 

• The annihilation diagrams of the four B — > a^K channels are similar with each other, the dominant contributions 
are all from the S — > P time-like form factor mediated by an (S + P) density. 

• B a — > (IqK and B~ — > a^if ~ are more complicated due to the appearance of the ao-emission diagrams. Because 
of the vanishing vector decay constant of ao, the factorizable emission diagrams are zero. For the nonfactorizable 
diagrams, the QCD penguin operators cancel for the neutral state of isospin triplet. The electroweak penguin 
operators have a small Wilson coefficients, thus the emission contributions in these channels are small. For 
the tree operators, although they are suppressed by the CKM matrix elements, the nonfactorizable emission 
diagrams can be enhanced for the large Wilson coefficients C2. So it is expected a large CP asymmetry in the 
decays B° -> a°if° and B~ -> a%K~ . 

From the above discussion, we can see that the dominant contributions are from the annihilation diagrams and the 
diagrams with tree operators. 



IV. RESULTS AND DISCUSSIONS 



For numerical calculations, we have employed the parameters in Tab. I. For the B meson wave function, we adopt 
the Gaussian- type model [14] (we choose the shape parameter oj — 0.4). As for the light-cone distribution amplitudes 
(LCDAs) of the pion and kaon, we use the results from QCD sum rules up to twist-3 [24] . Other parameters relevant 
to the scalar mesons have been given in the second section. 
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TABLE I: Input parameters used in the numerical calculation 



Masses 


niK' = 1.412 GeV, 
m ao = 0.98 GeV 


rriQ —1.7 GeV, 
M B = 5.28 GeV 


ml 


= 1.4 GeV 


Decay constants 


f B = 0.19 GeV 


f K = 0.16 GeV 


h = 


0.132 GeV 


Life Times 


t b ± = 1.671 x 10~ 12 s 


t b o = 1.536 x 10~ 12 s 






CKM 


V tb = 0.9997 
V ua = 0.2196 


Vts = -0.04, 
V ub = 0.00367e- i60 ° 







TABLE II: Various decay amplitudes (xlO _2 GeV 3 ) in decay B~ -> K* n~ and B~ -> 





F e L (a 4 ) 


F e > 6 ) 


A^(C 3 )+M«(C B ) 


Ftifu) + Fa(a 6 ) + M L a {C s ) + M*(C B ) 




scenario I 


0.98 


-12 


0.78 + 1.2i 


7.4 + 13.3i 


8.8 - 11. Si 


scenario II 


-1.4 


17.9 


2.1 - 0.37i 


-5.8 - 17.8i 


-13.6 - O.llz 


B" a" X 


9.3 


-11.3 


0.09 - 0.74i 


2.0-9.0i 


8.6 + 1.2i 



A. The Branching Ratios and The CP Asymmetries 

Using the parameters in the above, we give the numerical results for different amplitudes of B~ — > K^n^ in Table 
II. The numerical results confirm that the emission diagram of (V — A)(V — A) operators Fg^ v (a,4) indeed give small 
contributions because of the small vector current decay constant. The (V — A)(V + A) operators give the dominant 
contribution F^{a§). The opposite sign between scenario I and scenario II comes from the decay constant of the K$ 
meson. The non-factorizable contribution M%^(Cz) is small due to the small Wilson coefficient. .M^^^C^) is 
even smaller because of the cancellation between the two diagrams. 

According to PQCD power counting, the annihilation diagrams are power suppressed, but the suppression is not 
so effective in some cases, such as when chiral enhancement existing. Usually there is a large imaginary part in the 
amplitudes of the annihilation diagrams, which is the source of strong phase in PQCD approach. The numerical 
results in Table II indicate that the annihilation diagrams in scenario I are more important than that in scenario 
II. There are also tree operators contributing to the annihilation diagrams F^ (02) + A4^(C\), which are Cabibbo 
suppressed, but they are essential in direct CP violation. 

Now it is straightforward to obtain the results for the CP-averaged branching ratio of B~ — > Kq°tt~ , which is 
given in Tab. III. Comparing the two scenarios, we find: the results from scenario II are twice as scenario I, the most 
important reason is the larger scalar decay constant in scenario II; secondly, the nonfactorizablc diagrams are small, 
they only change the branching ratio slightly; the annihilation diagrams play an important role in both scenarios, it 

TABLE III: Branching ratios for the decays B — » Kqw and B aoK(in units of 10~ 6 ). The first theoretical error is from the 
decay constant of the scalar meson, the second and the third one is Gengebauer moments Bi and Bs, the uncertainty caused 
the CKM angle 7 is very small which is not listed here. The experimental data listed here are the world average values by the 
Heavy Flavor Averaging Group (HFAG) [11]. 



Channel 


scenario I 


scenario II 


exp. 


Channel 


scenario I 


exp. 


B~ -> Kfn- 
B° -> K*-*+ 
B° -> K~f\° 
B- K*-tt" 


on 7+4.3+0.8+1.8 
zu - ' -3.9-0.8-1.6 

9n n +4.2+0.8+1.6 
zu - u -3.8-0.7-1.5 
ln n +2. 1+0.4+1.0 
lu - u -l. 9-0.5-0. 9 
11 n+2.4+0.4+0.7 
1 - L -°-2. 1-0.3-0. 7 


A7 c+11.3+3.7+6.9 
-°- 10. 1-3.6-5.1 

ao n+iO- 2+3. 1+7.0 
^°- u -9.1-2.9-5.2 
10 4+4.4+1.5+4.0 
10 -^-3.9-1.4-2.9 
90 c+6.8+1.9+3.2 
zo -°-6.1-l. 9-3.5 


41.2 + 4.2 
46.6t^ 
25.5 + 9.9 


B~ — » K° a - 
B° -> K- a + 
B° - K°a° 
B~ -» K-al 


a q+0. 8+1. 1+2.0 
u - 3 -0.7-l.l-1.7 
n 7+1.1 + 1.6+2.7 
J - ' -1.0-1.4-2.2 
a 7+0.5+0.7+1.1 
' -0.5-0.8-1.1 
q c+0. 4+0. 4+1.0 
°-°-0. 4-0.6-1.0 


< 3.9 

< 1.6 

< 7.8 

< 2.5 
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TABLE IV: Ratios of the branching fractions in B — ► Kqtt. For the experimental values, we only use the central values. 



Ratios 


isospin limit 


Scenario II 


Scenario I 


Experiment 


Ri 


0.5 


0.43 


0.48 


0.39 


R 2 


0.5 


0.61 


0.55 




Ra 


1.0 


1.02 


0.95 


0.81 



can enhance the branching ratios about 50% in scenario I, and about 30% in scenario II. The current experimental 
data [11] is also listed in Tab. III. The large branching ratio is consistent with the results in scenario II, so scenario II 
is more preferable than scenario I, this conclusion is consistent with [12]. But the difference with [12] is: we directly 
calculate the annihilation contribution in B — > Kq-k, rather than fit the B — > Kq-it data and then use the B — ► clqK 
data to rule out scenario I. 

The decay amplitudes for the B~ — ► K° decay are also listed in table II, the results indicate that the emission 
diagrams almost cancelled out, as expected in section III. The branching ratio is dominated by the annihilation 
diagrams which is at the same level as the B~ — > K^tt~ , and the induced branching ratio is about twice larger than 
the experimental upper bound. The results also show that scenario I is not supported by the current experimental 
data. 

In the above discussion, we concentrate on B~ — ► Kq°tt~ and B — ► K° . The dominant contribution in B — > Kqtt 
for other channels is the same with B~ — ► Kq°it~ by isospin symmetry, except the contribution from the electro-weak 
penguin and the tree operators which can violate isospin symmetry. To explore the deviation of the isospin limit, it 
is convenient to define the parameters below: 

_ B(B ^K^tt ) 

1 ~ B{B0^K*-w+Y 

B(B- -> K*-n°) 

2 B(B- - K^ir-)' 

_ r(B°) B(B~ -+ KS°ir-) 

3 r(B-) B(B0 ^ K^-n+Y { ' 

For B — > a K , the definition is the similar except Kq — > K, ir — > a . These parameters are the ratios of the branching 
fractions, which should be less sensitive to many nonperturbative inputs than the branching fractions, thus it is more 
persuadable to test these parameters. In the isospin limit, if we ignore the CKM suppressed tree diagrams and 
electro- weak penguins, R\, R2 and R3 should be equal to 0.5, 0.5 and 1.0. The deviations reflect the magnitude of 
the tree operators and the electro-weak penguins directly. Our results and the experimental data are given in table 
IV, where we use the central values in table III for the experiment data. In both scenarios, the deviations from 
isospin limit are not large, which shows that the QCD penguin are dominant in the branching ratios, both in emission 
diagrams and the annihilation diagrams. The three ratios for B — > a^K decays arc: R\ = 0.48, Ri = 0.51, R3 = 0.68. 
There is a large deviation for the ratio R3, and the reason is the large tree contribution. So the large direct CP 
asymmetry for B° — > K + is also expected. 

The contribution from different effective operators shown in Eqs. (31-41) have been categorized to two groups 
according to the different CKM matrix elements: 

A=^[V ub V; s T-V tb V*P], (49) 

where T/P denotes the amplitude which comes from the tree/penguin operators respectively. The charge conjugate 
channel decay amplitude is the same as Eq. (49) except the sign of the weak phase. The formula for the direct CP 
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TABLE V: Direct CP asymmetries (in units of %) 



Channel 


Scenario I 


Scenario II 


exp. 


channels 


A C p 


b~ ->Tq°w- 


-1.5 


-1.7 


°-8 


B~ -> ^°a" 


4 


B° - K;-n+ 


9.2 


0.22 


-7± 14 


S° -> ^"0+ 


-70 


B° - i?f tt 


-9.0 


-6.8 


-34 ± 19 


S° - K°ag 


-17 


B" K*-7r° 


16.0 


3.5 






-70 



asymmetry reads: 



= l^l 2 - |A| 2 = 2z sin 7 sin J 
CP |i| 2 -|A| 2 l-2zcos 7 cos(5 + z 2 ' 1 ' 



where z = \ y tb ^° ||y |, <$ is the relative strong phase between two groups of contributions. This equation indicates 
that the direct CP violation depends on the ratio of the tree and penguin contribution. The direct CP asymmetry 
is very small if the ratio is too large or too small, while the comparable tree and penguin contributions imply large 
direct CP asymmetries. For B — > K^tt, the penguin operators give the dominant contribution, but the tree operators 
suffer from the CKM suppression, so it is expected the direct CP asymmetry is small. We list our results in table 
V as well as the experimental data [11], where the results are consistent with the experiments in both scenarios. But 
in the decay B — > clqK, as mentioned above, the emission penguin contribution cancels, while the tree operators are 
large, sizable direct CP asymmetries are predicted in these channels, especially in B° — > K~ and B~ — > a°K~ . In 
QCDF, the central values of direct CP asymmetries for all four channels are very small, but with large uncertainties 
for B~ — > clqK~ and B° — > a^K~ . Furthermore, we may expect the similar size of CP asymmetries in similar decays 
B -> 00(1450)^. 



B. The Theoretical Uncertainties 



In the above calculation, the uncertainties from the decay constants can give sizable effects on the branching ratio, 
but not to the direct CP asymmetries. Furthermore, there are other sources of uncertainties: 

• The twist-3 distribution amplitudes of the scalar mesons are taken as the asymptotic form for lack of more 
reasonable results from non-perturbative methods. These distribution amplitudes will be studied in the future 
work [25]. In [25], we find the Gegenbauer moments of the twist-3 distribution amplitudes are rather small, 
which implies the results will not be changed sizably. 

• The Gegenbauer moments B\ and Bz for twist-2 LCDAs of Kq and ao have sizable uncertainties, which can 
lead to the theoretical errors. We include these uncertainties in the results and they can give about 20% <~ 30% 
uncertainties to the branching ratio. 

• The uncertainties of the light pseudoscalar meson and B meson wave functions, the factorization scale, et al. 
have been studied extensively in [26]. The uncertainty from the factorization scale is within 10%. The major 
source of the uncertainty comes from the meson distribution amplitudes. The results can be varied by (10 — 30)% 
by changing the parameter in the wave functions. 

• The sub-leading order contributions in PQCD approach have also been neglected in this calculation, but these 
corrections have been calculated in refs. [27] for B —> irK, irir, etc. These corrections can change the penguin 
dominated processes about 20% of the branching ratio of B — > irK. We may expect similar effect in B — > KqIt. 
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• The uncertainties of CKM matrix elements and the CKM phase angle can also affect the branching ratios and 
CP asymmetries. In the two kinds of decays considered in this paper, the decay amplitudes are the functions of 
the CKM angle 7, whose value given in PDG06 is 7 = (63t^)° [28]. With the CKM angle 7 varying at this 
area, we find that the error area is very small in B — > K^tt decays in both scenarios. For the B — ► clqK decays, 
the error area is some larger, but within ten percent. 

• The long distance re-scattering can also affect the branching ratios and CP asymmetries. This effect could be 
phcnomenologically studied in the final-state interactions [29]. We need more data to determine whether it is 
necessary to include the re-scattering effect in B — > SP decays. 

V. SUMMARY 

In this paper, we calculate the decay modes B — > KqTt and B — > ciqK within pcrturbativc QCD framework. For 
B — > Kq-jt, we perform our calculation in two scenarios of the scalar meson spectrum, our calculation indicates that: 
scenario II is more consistent with the experimental data than scenario I. We directly calculate the contribution from 
annihilation diagrams: it can enhance the branching ratios about 50% in scenario I, and about 30% in scenario II. 
Our predicted branching ratio of B — > a K in scenario I is larger than the experimental upper bound, which indicates 
ao(980) can not be interpreted as qq. We calculated the direct CP asymmetries and the isospin parameters in these 
decays, and we find that in B — > K^ir (in both scenarios) the direct CP asymmetries are small, which are consistent 
with the present experiments; the deviation from the isospin limit is also small. There is large CP asymmetries in 
B — > a n (980)K due to the relatively large tree contributions in scenario I. We expect similar CP asymmetries in 
B -» 00(1450)^ . 
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